Cluster model of decagonal tilings 
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A relaxed version of Gummelt's covering rules for the aperiodic decagon is considered, which 
produces certain random-tiling-type structures. These structures are precisely characterized, along 
with their relationships to various other random tiling ensembles. The relaxed covering rule has 
a natural realization in terms of a vertex cluster in the Penrose pentagon tiling. Using Monte 
Carlo simulations, it is shown that the structures obtained by maximizing the density of this cluster 
are the same as those produced by the corresponding covering rules. The entropy density of the 
covering ensemble is determined using the entropic sampling algorithm. If the model is extended 
by an additional coupling between neighboring clusters, perfectly ordered structures are obtained, 
like those produced by Gummelt's perfect covering rules. 
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I. INTRODUCTION 

Many quasicrystals are completely covered by overlap- 
ping copies of a single cluster (Fig. . Two overlapping 
clusters must agree in the overlap region, which restricts 
the possible relative positions and orientations of neigh- 
boring clusters. Cluster overlaps therefore create order, 
in favorable cases even perfect quasiperiodic order. 

This observation can be used to formulate several vari- 
ants of an ordering principle for quasicrystals (for a re- 
view, see Ref. 0): A perfect quasicrystal can be obtained 
by requiring either that a given cluster completely cov- 
ers the structure, or that the cluster has maximal den- 
sity in the structure, or that it covers the structure with 
maximal density. With such ordering principles, perfect 
quasiperiodic order could be obtained for decagonal, 4,5 
octagonal^ and dodecagonal£ tilings and quasicrystals. 
Assuming that such a cluster is an energetically preferred 
atomic configuration, the maximization of the cluster 
density minimizes the free energy. With this hypothe- 
sis, the covering approach might serve as a simple ther- 
modynamic mechanism for the formation of quasicrys- 
tals. The covering approach can be regarded as a partic- 
ularly simple realization of energy based matching rules, 
where only the most important local configurations (the 
clusters) have to be preferred energetically^, not all al- 
lowed local configurations. As a variant, it has also been 
suggested 9 to penalize the worst local configurations, in- 
stead of prefering the best ones, which provides another 
way to simplify the matching rule approach. 

The same ordering principles can also be used to pro- 
duce supertile random tiling structures, which are locally 
ordered but show disorder on larger scales. These struc- 
tures are obtained whenever the chosen cluster is not 
selective enough and hence allows too many different 
overlaps^ This happens, in particular, if the cluster is 
too small to restrict the number of different overlaps ^ 
or if it is too symmetric. In this respect, it is interesting 
to note that asymmetric clusters seem to be preferred 



by the electronic structure in decagonal quasicrystals»ii 
Whereas supertile random tilings are locally ordered for 
energetic reasons, their long range order is produced by 
entropy maximization, as is the case for other random 
tilings 1 ^. 

In this paper, we will concentrate on cluster models 
for decagonal structures. Gummelt's aperiodic decagon 4 
provides a striking example how perfect quasiperiodic or- 
der can be obtained by a simple cluster covering princi- 
ple. This example has been so convincing, that many 
researchers tried to map their experimental structures to 
the Gummelt decagon, even though the fit in the overlap 
was often not perfect (see, e. g., Refs . ITTI and fl^) . and the 
overlapping constraints not exactly equivalent. However, 




FIG. 1: High-resolution transmission electron microscopy im- 
age of decagonal Al-Ni-Co (courtesy of S. Ritsch and C. Beeli, 
compare Ref. U|). The superimposed tiling has been recon- 
structed by an automated procedure^ 
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many experimental decagonal quasicrystal structures are 
not perfectly quasiperiodic, and it is therefore interest- 
ing to consider also overlap rules which are less restric- 
tive than the perfect rules of Gummelt, and which do not 
enforce perfectly ordered, but rather (supertile) random 
tiling structures. 

The analysis of such relaxed overlap rules and their 
corresponding structures will be the main topic of this 
paper Mi In Sec.|n] two different relaxed versions of Gum- 
melt's overlap rules are discussed, and the structures 
which they produce are precisely characterized, along 
with their relationships to various other random tiling 
ensembles. Subsequently, in Sec. IIIII we introduce a ver- 
tex cluster in the Penrose pentagon tiling (PPT) whose 
structure imposes the previously discussed overlap con- 
straints in a natural way. It is shown by Monte Carlo 
(MC) simulations that the structures with maximal den- 
sity of this cluster are the same as those produced by the 
corresponding overlap rules. In Sec. IIVI we determine 
the entropy density of the set of states with maximal 
cluster density, using the entropic sampling algorithm. 
An additional coupling between neighboring clusters is 
introduced in Sec. [3 and it is shown that this coupling 
is capable of ordering the random tilings to perfectly or- 
dered structures. 



II. COVERINGS FOR PERFECT AND 
RANDOM PENROSE PENTAGON TILINGS 

It is well known that each covering of the plane by 
Gummelt's aperiodic decagon (Fig. |2a)) is equivalent 
to a perfect Penrose tiling, 4 if the covering has the fol- 
lowing property: Whenever two decagons overlap, their 
colorings agree in the entire overlap region. It turns out 
that Gummelt's rule permits only two different types of 
overlaps, which are shown in Fig. |2b): the smaller A- 
and the larger B-overlaps. Furthermore, due to the col- 
oring, there are only certain overlap zones for allowed 




FIG. 2: Gummelt decagon (a), representative A- and B- 
overlap (b), and allowed overlap zones (c). 




FIG. 3: Perfect PPT, superimposed on the corresponding 
Gummelt decagon covering. 



overlaps with neighboring decagons: four for A- and two 
for B-overlaps (Fig.EIc)). This altogether is what we will 
call the perfect rule (in order to distinguish it from other 
variants being discussed later). 

The decagon centers of such a perfect covering form 
the vertex set of a perfect PPT (Fig.|3). Conversely, each 
PPT can be obtained from exactly one covering satisfying 
the Gummelt overlap rules. We therefore have a local 
one-to-one correspondence between PPTs and Gummelt 
coverings. As the Gummelt decagon represents a cluster 
in the corresponding quasicrystal, we will often use the 
term cluster for the covering decagon. 

In order to allow for partially disordered coverings, 
Gummelt et alM^ have proposed to relax the over- 
lap rules to some extent. To understand the type of 
relaxation, recall that if the perfect rules are obeyed, 
a decagon may have small A-overlaps with neighbor- 
ing decagons in four possible directions, and bigger B- 
overlaps with neighboring decagons in two possible di- 
rections (Fig. |2Ic)). The coloring in the overlap region 
has an orientation, which must be respected. All possi- 
ble overlaps are therefore oriented. As a relaxation of the 
perfect rule, Gummelt et ali 15 i 16 have proposed to aban- 
don this orientation constraint, and to retain only the 
non-oriented overlap zones, as shown in Fig. |^c). This 
overlap rule will be referred to as the fully relaxed rule. 

There is a natural intermediate rule between the per- 
fect and the fully relaxed rule. In this variant, which 
will be called the relaxed rule, the orientation condition 
is abandoned only for the small A-overlaps, but is re- 
tained for the larger B-overlaps. This kind of overlap 
rule can be motivated physically as follows: The large 
B-overlaps result in a strong interaction between the two 
overlapping clusters, which must be in its ground state, 
whereas the small A-overlaps only lead to a weak inter- 
action with a small energy difference between differently 
oriented A-overlaps. This intermediate rule and the re- 
sulting structures will be the main topic of this paper. 
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FIG. 4: Random PPT with all the spiky tiles (shaded in gray) 
surrounded by pentagons. This tiling is equivalent to a ran- 
dom HBS tiling (gray lines), whose tile edges connect the 
centers of neighboring pentagons. 

Gummelt et ali^i^ have shown that each covering sat- 
isfying the fully relaxed rule has the property that its 
cluster centers form the vertex set of a random PPT. It 
has the additional property that all the spiky tiles (stars, 
ships, and rhombi; shaded in gray in Fig. are com- 
pletely surrounded by pentagons. (In the following, when 
we use the term "random PPT" , we always mean one 
satisfying this extra condition; more general ones do not 
play any role here.) Such a random PPT is equivalent to 
a random hexagon-boat-star (HBS) tiling (gray lines in 
Fig-EJ. Since coverings satisfying the more restrictive re- 
laxed rule also satisfy the fully relaxed rule, their cluster 
centers form the vertex set of a random PPT, too. Con- 
versely, it is easy to see that every random PPT can arise 
both from relaxed and from fully relaxed coverings. The 
only difference between relaxed and fully relaxed cover- 
ings is the number of coverings associated with a given 
random PPT. 




FIG. 5: The orientation of a decagon on a vertex in a PPT 
is fixed in the case of four A- or two B-neighbors (left). The 
only vertices with a choice for the decagon orientation are the 
obtuse rhombus corners, which have three A-neighbors and 
one B-neighbor (right). 



To see this, we note that the orientation of a cluster on 
a vertex in a PPT is completely fixed by the presence of 
four A-neighbors or two B-neighbors (Figs.[S]and[fJa-d)), 
as in these cases the four A- or the two B-overlap zones, 
respectively, are completely saturated. A-neighbors are 
separated by an edge of a tile or a long diagonal across 
a ship or star, whereas B-neighbors are separated by a 
short diagonal of a rhombus, ship or star. The only ver- 
tices whose cluster orientation is not fixed by their lo- 
cal environment in the tiling are the obtuse corners of 
the rhombi, which have three A-neighbors and one B- 
neighbor (Figs.[5]and[Ue)). Therefore, two cluster orien- 
tations are possible for each obtuse rhombus corner. For 
the fully relaxed rule, where no orientation conditions for 
the overlaps have to be obeyed, we thus have altogether 
four choices per rhombus for the cluster orientations. 

However, for the relaxed overlap rule we have to obey 
the orientation condition for the B-overlaps. It is eas- 
ily shown that this condition is always satisfied for B- 
overlaps across ships and stars. In order to fulfill the 
orientation condition for B-overlaps across rhombi, too, 
the orientations of the clusters on opposite obtuse rhom- 
bus corners cannot be chosen independently. If for one 
of the corners an orientation is chosen, the orientation of 
the other is already fixed, i. e., the condition is satisfied 
only for two of the four possible combinations mentioned 
above. 




FIG. 6: Orientation of a decagon on a vertex in the PPT, 
depending on the local environment. The orientations are 
given by an arrow, as defined in the upper right corner. 
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FIG. 7: Relationship between 4-level (thin black lines) and 
2-level random tilings (thick gray lines). The matching rules 
for the "double" arrows (here drawn in black) are still obeyed. 
Since the matching rules for the "single" arrows (here drawn 
in white) are no longer maintained, each hexagon of the 2- 
level or random HBS tiling can be subdivided in two different 
ways, as shown in the bottom right corner. 

In the same way, one can quantify the relationship be- 
tween the cluster coverings and certain variants of ran- 
dom Penrose rhombus tilings. The random HBS tilings 
arise from random Penrose rhombus tilings still satisfying 
the matching rules for "double" arrows (drawn in black 
in Fig. UJ . Such random Penrose rhombus tilings are also 
called 4-level random tilings. 12,17 When the edges with a 
double arrow are simply wiped out, we obtain the random 
HBS tilings, which are also known as 2-levcl tilings. 12 The 
relationship between 4-level and 2-level random tilings is 
not one-to-one: Whereas the subdivision of boats and 
stars is unique, there are two choices for the subdivision 
of each hexagon into rhombi (Fig. 0), just as there are 
two possible cluster assignments on the obtuse rhombus 
corners in the PPT, as discussed above (Fig.HJe)). Since 
rhombi in the PPT and hexagons in the HBS tiling are in 
one-to-one correspondence (Fig. @J, this implies that the 
multiplicity of relaxed cluster coverings and 4-levcl ran- 
dom tilings, related to a given random PPT, is the same. 
Apart from the extra multiplicities, the relaxed cover- 
ing rule is therefore equivalent to the Penrose double- 
arrow matching rules (ignoring single arrows), whereas 
the Gummelt convering rule is equivalent to the full Pen- 
rose matching rules. 

III. CLUSTER DENSITY MAXIMIZATION 

In the last section, we have considered cluster cover- 
ings, where our clusters have simply been decagons with 



certain overlap rules. Another variant of an ordering 
principle for quasicrystals based on a cluster picture is 
the cluster density maximization^^ which we will con- 
sider in the following. 

The relaxed overlap rule discussed above allows for a 
very natural realization in terms of a vertex cluster in the 
PPT. This vertex cluster is shown in Fig. |S| We have to 
point out that the tile edges are drawn only as a guide to 
the eye; they are not part of the cluster, only the vertex 
set counts. It is easy to sec that the vertex set of the clus- 
ter cannot enforce the orientation of the small A-overlaps 
of the Gummelt decagon, whereas the orientation of the 
B-overlaps is intrinsically enforced. The A-overlap con- 
sists of a rhombus or a rhombic area inside a ship or star 
without orientation, and the B-overlap is formed essen- 
tially of a hexagon-shaped area with an interior vertex in 
asymmetric position which yields an orientation (Fig.|SJ). 

With this cluster, we can build a statistical model for 
the cluster density maximization. We consider the set 
of all random PPTs (we still require that spiky tiles are 
completely surrounded by pentagons) and assign to each 
tiling a statistical weight which is simply the number 
of vertex clusters it contains. With a suitable MC algo- 
rithm, it is then possible to find the subensemble of those 
random PPTs which have maximal cluster density. For 
this purpose, we need a MC dynamics which is ergodic 
in the ensemble of all random PPTs. By repeated flips it 
is then possible to turn any random PPT into any other. 

We have found (see Sec. HV|) that the flip moves shown 
in Fig.[5]have the required properties. The flip configura- 
tions consist of a hexagon with an interior vertex which 
can jump to its "mirror image" . This move corresponds 
to a change in the orientation of the hexagon by 180°. 
The vertices of the hexagon itself are not affected, but 
the adjacent tile configurations are changed depending 
on the local environment. In Fig. EJa), e. g., the adja- 
cent ship and rhombus are exchanged, or in Fig. E£b), 
a star and a rhombus are transformed into two ships, 
etc. These flips preserve the property that the spiky tiles 
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FIG. 8: Vertex cluster, superimposed on the Gummelt 
decagon (top left), and representative A- and B-overlaps (top 
middle and right). This cluster enforces the relaxed overlap 
rules. The bottom row shows examples of tile configurations 
for the vertex cluster. 
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FIG. 9: Flip moves for the MC simulation. The flips (a- 
c) are allowed, whereas move (d) is forbidden since it would 
produce a new kind of tile (the zigzag-shaped configuration 
at the bottom). 



are always surrounded by pentagons. There are some 
"flip configurations" where the local environment pro- 
hibits the flip. This is the case for configurations like in 
Fig. Eld), where the flip would introduce a new kind of 
tile (see also Ref . IT^l . In all our simulations, such flips 
were forbidden. For the sake of completeness, it is explic- 
itly shown in Fig. 1101 how a new cluster can be created 
by a single flip. 

The flips used for the PPT are in direct correspon- 
dence with hexagon flips in the associated 4-level Pen- 
rose tiling. This is illustrated in Fig. ^1 Note that the 
other type of hexagon flip plays no role, because it leaves 
the HBS tiling, and thus the PPT invariant. As it is 
known that hexagon flips are ergodic for the 4-level Pen- 
rose tiling, this correspondence adds further confidence 
that the PPT flips are indeed ergodic, too. 

With this MC scheme, the states of maximal cluster- 
density can be determined by simulated annealing, using 
as energy the negative of the number of clusters, thus 
mimicking the total cohesion energy of the clusters. Then 
the energetic ground state, reached at low temperatures, 
is the ensemble of states with maximal cluster density. 
The method we use is based on the Metropolis impor- 




AE = -l 



FIG. 10: Creation of a new cluster by a single flip, which 
lowers the total energy. 






(c) 






FIG. 11: Relationship between the flip configurations in the 
PPT (black), as shown in Fig. [5] and the corresponding lo- 
cal tile configurations in the random Penrose rhombus tiling 
(gray), which were used in Ref.ll7l The flip move is explicitly 
shown in the bottom right corner. 



tance sampling algorithmiiLiS The basic MC move is as 
follows: (i) Choose a vertex randomly, (ii) If it can be 
flipped, calculate the energy change AE (which is the 
negative of the change in the number of clusters) and flip 
it with probability 



e -PAE for AS > , 

1 otherwise . 



(1) 



This algorithm fulfills the condition of detailed balance. 

It turns out that the states of maximal cluster den- 
sity are supertile random PPTs, whose tiles have an 
edge length t 2 times that of the small tiles (where 
t = (1 + \/E)/2 is the golden number). An example 
of such a supertiling is shown in Fig. 1121 It cannot be a 
perfect tiling, since it is still possible to move clusters in 
the ground state without changing their number. This is 
due to the fact that the relaxed rule and hence our vertex 
cluster does not enforce the orientation of the A-overlaps. 

In view of the results of the previous section, this is 
of course not too surprising. The cluster centers sit on 
the vertices of the supertiling, covering all vertices of the 
small tiles. Since the vertex cluster is smaller than the 
Gummelt decagon, it does not cover the whole area, but 
only the vertices; there remain small pentagons uncov- 
ered, which sit at the center of the supertile pentagons. 
This does not affect the overlap constraints, however. 



FIG. 13: In the supertile random PPT, each rhombus is 
counted twice because of the two possible cluster orientations 
on the obtuse corners (left). This is in one-to-one correspon- 
dence with the hexagons in the 4-level random tiling (right). 



FIG. 12: Structure with maximal cluster density. The cluster 
centers form the vertices of a supertile random PPT (gray 
lines). 



Our results therefore imply that there is a one-to- 
one correspondence between decagon coverings satisfy- 
ing the relaxed rule and structures with maximal den- 
sity of the vertex cluster (Figs. 151 and 1121 respectively). 
Although these two ordering principles are very similar, 
they are conceptually slightly different and have to be 
distinguished. 

Supertile random tiling ensembles as a result of cluster 
maximization were found already in Ref. 0. Maximizing 
star decagon clusters in a rhombus tiling leads to an en- 
semble of HBS-type supertile tilings*^ That ensemble 
contains also other structures, however, which is not the 
case for our cluster. In particular, it should be noted that 
if the tile stoichiometry of the pentagons and spiky tiles 
admits an HBS supertile tiling, then the state of maximal 
cluster density is always an HBS supertile tiling. A phase 
separation as discussed in Ref. is not possible. The only 
tiles that could be separated are the thin rhombi, but if 
there are enough pentagons, it is always advantageous 
to surround the rhombi with pentagons. The ensemble 
obtained by maximizing our vertex cluster is therefore 
a strict subensemble of the one obtained by maximizing 
the star decagon in random rhombus tilings. 9 



IV. ENTROPY DENSITY 

With our cluster model, it is also possible to measure 
the entropy density of the ensemble of structures with 
maximal cluster density and thus the entropy density of 
the relaxed cluster covering ensemble. In the previous 
section, we have introduced an energy model which as- 
signs a cohesion energy to each cluster in the structure. 
In this model, the ground state, i. e., the state of max- 
imal cluster density, consists of supertile random PPTs 
with an extra weight of two per rhombus, because for 
each rhombus there are two choices of a cluster config- 
uration with the same number of clusters, as shown in 



Fig- El ( see also Fig. E{e)). At infinite temperature, on 
the other hand, we have the full random PPT (at the 
level of the small tiles), where each rhombus is counted 
only once. 

With the entropic sampling algorithmji2i2£ we can de- 
termine the entropy of the system as a function of energy. 
In this method, the Boltzmann probability e~@ E of the 
Metropolis algorithm is replaced by the factor e~ s ( E \ 
where S(E) is the microcanonical entropy function. Since 
this factor is just the inverse of the number g(E) of states 
with energy E, according to 



S(E)=\ng(E) 



(2) 



(in units of /cb), we obtain a uniform energy distribution 
P{E) oc g{E)e- s(E) = 1 , (3) 

which corresponds to a random walk through the energy 
space of the system. 

However, the exact entropy function S(E) is not known 
a priori, hence it has to be determined iteratively, start- 
ing with a rough estimate of S(E). This estimate can 
be obtained by a short run using as "entropy function" 
S(E) = 0, which yields an estimate for the degeneracy 
g{E) of the different energies and thus, via Eq. an es_ 
timate for the real entropy function. Another possibility 
to get a good estimate is to take advantage of the exten- 
sive nature of entropy by scaling the entropy function of 
a smaller system to a larger one. 

Subsequently, the entropy function is optimized by an 
iterative procedure. Analogously to the Metropolis algo- 
rithm, we choose for the flip probabilities 



e~ A5 forA5>0, 
1 otherwise , 



(4) 



where AS" = S(E + AE) — S(E) is the change in entropy 
due to the considered MC move (with energy change 
AE). This choice of p likewise fulfills the condition of 
detailed balance. The iteration scheme is then as follows: 
(i) Run a MC simulation based on Eq. |0J in order to ob- 
tain an energy histogram H(E). (ii) Use the measured 
histogram to correct the entropy function according to 
the update rule 



S(E) 



S{E) +lnH(E) for ^0, 
S(E) for H(E) = . 



(5) 
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FIG. 14: Representative example of the microcanonical en- 
tropy as a function of the energy of the system. 

(iii) Continue at (i) until the energy histogram is suffi- 
ciently uniform. 

An example of such an entropy function is shown in 
Fig. [21 As this method does not yield absolute entropy 
values, we can only measure entropy differences, in par- 
ticular the difference S between the ground state and 
the infinite-temperature state (i. e., the maximum of the 
entropy function). The entropies at zero (So) and in- 
finite temperature (5*00) are both entropies of random 
PPTs, once with a twofold degeneracy for each rhombus 
(Fig. I13[) . and once without. To compare these two en- 
tropy values, the latter one has to be corrected by adding 
an extra double-counting of the rhombi, which yields for 
each rhombus a factor of 2 in the degeneracy g^ or an 
additive contribution of In 2 in the entropy , respec- 
tively. Moreover, the two random tilings are at different 
length scales, since the supertile edges are t 2 times larger 
than the small tile edges, which has to be taken into ac- 
count due to the extensive nature of entropy. Therefore, 
we have the following relation between the two entropy 
densities ao and a^: 

T 4 a = doc + p T h In 2 , (6) 

where p r h is the measured rhombus density in the infinite- 
temperature state. If we write = ao + a, we end up 
with an equation for the ground state entropy density ao, 
in which all other quantities can be measured: 

go = 4 1 : (o- + prh In 2) . (7) 

T 4 — 1 

The ground state entropy density has been determined 
in this way for several periodic approximants. By finite- 
size scaling, the values can then be extrapolated to infi- 
nite system size. For this purpose, the entropy density, as 
a function of system size N, is fitted to a function of the 
(empirical) form a + bc~ c / N , which proves to work very 
well fFig. H5fl . At a scale where the supertile edges (which 
separate A-overlap neighbors in the cluster model) have 
unit length, we obtain a value of 

a /k B = 0.253 ±0.001 (8) 
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FIG. 15: Finite-size scaling of the entropy density of the re- 
laxed cluster coverings, where N is the number of vertices. 
The applied fit function is of the (empirical) form u(N) = 
a + be~ c/N . In this graph, the line log[cr(jV) - a] = \ogb-c/N 
is shown. The error bars are smaller than the plot symbols. 

for the entropy density of the relaxed coverings. This 
can be compared with the value which Tang and Jaric 
have obtained by Metropolis-type MC simulations for the 
entropy density of the 4-level random tiling^ In Sec.lTTl 
we have seen that 4-level random tilings are in one-to-one 
correspondence with relaxed cluster coverings (see also 
Fig. 1131) . If the different length scales of the two tilings 
are taken into account by a simple geometric conversion, 
the value of Tang and Jaric turns into 

a /k B = 0.255 ±0.001, (9) 

which is compatible with our result. 

V. COUPLING BETWEEN CLUSTERS 

The only difference between the perfect and the relaxed 
overlap rule is that the latter does not require oriented 
A-overlaps (whereas the orientation of the B-overlaps has 
to be obeyed). Since not all relaxed coverings are per- 
fect, there must be A-overlaps which do not obey the 
orientation condition of the perfect rule. A closer anal- 
ysis showsi£*i& that there is actually only one kind of 
disoriented A-overlap. All A-overlaps which can occur 
in relaxed coverings or supertile random PPTs, respec- 
tively, are shown in Fig. El For the disoriented A-overlap 
not permitted by the perfect rule (Fig. UST dll. the two 
clusters have antiparallel orientations. 

To order the (supertile) random tiling structures to 
perfect tilings, we introduce a coupling between neigh- 
boring clusters in such a way that overlaps which are not 
permitted by the perfect rule are energetically penalized. 
We expect such a coupling to be weak, because these 
kinds of defects can be detected only at larger scales. 
However, at low temperatures this coupling might still 
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FIG. 16: Possible A-overlaps in supertile random PPTs or 
relaxed coverings, (a-c) obey the perfect rule. Only (d) is a 
disoriented overlap. 




FIG. 17: Supertiling with cluster orientations indicated by 
arrows. The forbidden A-overlaps, corresponding to tile edges 
with antiparallel arrows at their ends, are marked. 



be able to order the supertile random tiling ground state 
of the relaxed cluster covering to a perfectly quasiperi- 
odic structure. 

This suggests a scenario with two energy or temper- 
ature scales. The presence of each vertex cluster lowers 
the cohesion energy by a large amount, so that structures 
with maximal cluster density are strongly favored, even 
at relatively high temperatures. The equilibrium struc- 
tures at these temperatures are therefore relaxed cluster 
coverings. Additionally, there is a small coupling between 




FIG. 18: Basic flips in the underlying tiling (top) and effective 
moves on the level of the supertiling (bottom). The hexagon 
flip (left) is the same as the one used before. Additionally, we 
also have the rhombus flip (right) which does not change the 
vertex structure of the supertiling. For both flip types, the 
cluster orientations have to be adjusted consistently with the 
underlying structure. 




FIG. 19: Creation (top) and shift of defects (bottom) by single 
flips of the two types: hexagon flip (left) and rhombus flip 
(right). 



neighboring clusters, which can order the supertile ran- 
dom tiling to a perfect tiling at low temperatures. 

We have verified the feasibility of this approach by 
MC simulations. For this purpose, we consider only the 
subensemble of states with maximal cluster density. In 
other words: We run our simulations at the level of the 
supertiling. If we represent a cluster in the supertile PPT 
by its center and its orientation (given by an arrow), we 
can describe the covering in a much more compact way, 
as shown in Fig. 1171 Such a setup keeps the number of 
clusters constant, so that we cannot leave the states of 
maximal cluster density, which simplifies the simulation 
considerably. 

As flip moves we can still use those of Fig. |51 except 
that we now have to adjust the cluster orientations of 
several vertices in the flip configuration. We have to do 
this in a way consistent with the vertex configuration in 
the underlying tiling, which actually determines the clus- 
ter orientations (Fig. ll8|) . The first type of flip consists of 
two simultaneous flips in the underlying tiling. We have 
to point out that for a hexagon-type flip configuration 
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like the one in Fig. [^J (bottom left) the cluster orienta- 
tions on the obtuse corners of the adjacent rhombus have 
to be as shown. Otherwise, the flip is not allowed, since 
then the cluster orientations would be inconsistent with 
the ship-shaped vertex configuration after the flip. 

Additionally, we have to introduce a new type of flip 
which only changes the cluster orientations on the obtuse 
corners of a rhombus, but keeps the tiling itself fixed. 
Such a move corresponds to a single flip in the underlying 
tiling. In comparison with the flip moves in a Penrose 
rhombus tiling, 12 ' 1 - the hexagon flip corresponds to D- 
type and the rhombus flip to Q-type configurations i^i 

This MC dynamics can change the number of defects, 
as demonstrated in Fig. Again, we have used a 

Metropolis-type MC schemefi2*i^ now with the number 
of defects as energy. We have seen in our simulations 
that the coupling of the clusters, which penalizes the de- 
fects, is indeed capable of ordering the random tilings to 
perfectly quasiperiodic structures. In other words: The 
ground state, reached by simulated annealing, is a perfect 
PPT, whereas the high-temperature state is a supertile 
random tiling or relaxed cluster covering. In this respect, 
"high" temperature means high compared to the cluster 
coupling, but still low compared to the energy required 
to break up clusters. 

We have to mention that the ground state is not al- 
ways a perfect PPT. Perfectness is determined by lifting 
the vertices to hyper-space; if the projections of the ver- 
tices onto the perpendicular space are all inside the ac- 
ceptance region, the tiling is perfect, otherwise it is not. 
Since we use periodic approximants of PPTs in our sim- 
ulations, the minimal number of defects is always larger 
than zero —2 These defects can be shifted without chang- 
ing the energy. In this process, a vertex can leave the 
acceptance region in the perpendicular space. Thus, the 
tilings with minimal number of defects are not necessarily 
perfect. 

With the model of coupled clusters, it is also possible 
to measure the entropy density of the relaxed covering en- 
semble. In this case, the ground state is ordered and has 
zero entropy (at least in the thermodynamic limit), and 
the high-temperature state is the one whose entropy we 
are interested in. We therefore only need to measure the 
difference between the entropies of the high-temperature 
state and the ground state, and extrapolate these values 
to infinite system size. It turns out, however, that the 
finite-size scaling for this model does not work as well 
as for the model with a supertile random tiling in the 
ground state as considered in Sec. IIVI The results are 
therefore less precise, but nevertheless consistent. 



VI. DISCUSSION AND CONCLUSION 

In this paper, we have discussed different ordering 
principles for quasicrystals based on the cluster picture, 
namely the cluster covering principle and the principle of 
cluster density maximization. A relaxed version of the 




FIG. 20: Atomic cluster found by Roth and Henley— in a 
molecular dynamics simulation. White atoms are at z = 0, 
black atoms at z = | , and dotted atoms at z = |, | (in units 
of the period in 2-direction) . The area of our vertex cluster is 
shaded in gray. At the center of a star (marked with a box), 
we would have expected a different atomic configuration (see 
text). 

covering rules for Gummelt's aperiodic decagon has been 
considered, which produces as ground state a supertile 
random PPT. It has been shown that this relaxed over- 
lap rule has a very natural realization in terms of a vertex 
cluster in the PPT. 

The feasibility of our model has been tested by MC 
simulations. In particular, we have verified that the re- 
laxed cluster coverings coincide with the states of maxi- 
mal cluster density. The entropy density measured in the 
random covering ensemble is found to agree with the en- 
tropy density obtained by Tang and Jaric for the equiva- 
lent random 4-level tiling. Moreover, we have shown that 
a coupling between neighboring clusters can order the 
random-tiling-type ground states to a perfectly quasiperi- 
odic structure. 

Models of this kind can be very suitable for the expla- 
nation of experimentally observed decagonal quasicrys- 
tals, as the latter are often not perfectly quasiperiodic, 
but more of a random tiling nature^ Random PPTs, 
and the closely related HBS tilings, are often observed in 
high-resolution electron microscopy images (Fig. ^ com- 
pare Figs. 151 and IT21) ±1 

The most striking resemblence, however, is with a 
three-dimensional atomic cluster found by Roth and 
Henley22i in a molecular dynamics simulation of binary 
decagonal Frank-Kasper-type quasicrystals (Fig. I20fl . 
whose lateral overlap constraints are almost the same 
as those of our two-dimensional vertex cluster. The only 
discrepancy is at the center of a star, where for full equiv- 
alence we would have expected a single atom at z = in 
the Roth- Henley cluster, not two atoms at z — |, |. 

The vertex set of our cluster can be realized by several 
tile configurations in the PPT (Fig. [SJ) • This might pos- 
sibly be related to the results of an Al-Co-Ni structure 
analysis by Cervellino, Haibach, and Steurer—4 They ob- 
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served that there is a perfectly quasiperiodic long-range 
order of the centers of atomic clusters, but the interior 
structure of the clusters themselves is disordered. This 
means that the local atomic interactions cannot enforce 
local order, but there exists a long-range order. The pro- 
posed mechanism for long-range ordering is the electronic 
long-range term of free electrons. These long-wavelength 
electrons "see" only a simpler, "smeared out" version of 
the complex cluster structure. 

As we have seen in Sec. [V] the shift of a flipping clus- 
ter center is large, whereas the corresponding two vertex 
moves in the underlying tiling are small (Fig. 118(1 . In 
the atomic simulations of Al-Co-Ni done by Honal and 
Wellberryji& the flip of a cluster corresponds to jumps of 
only four atoms. The magnitude of these jumps is small 
(about 1 A) compared to the scale of the clusters (about 
20 A). 

All our simulations have been in two dimensions. It is 
well known that, for hnite-range interactions, quasicrys- 
tal structures in two dimensions cannot be stable at any 
positive temperature. 25 At non-zero temperatures, they 



are always in a random tiling state, the "phase transi- 
tion" from the ordered phase to the random tiling state 
being at zero temperature. However, in three dimensions 
the critical temperature is expected to be positive. For 
this purpose, one can consider as three-dimensional sys- 
tem a stacking of our two-dimensional models with a suit- 
able coupling between neighboring layers. Such layered 
tiling models have already been studied by Jeong and 
Steinhardt but are possible also in the cluster cover- 
ing setting. As expected for three-dimensional systems, 
the phase transition from ordered to random-tiling-type 
structures is at finite temperature. The results of these 
studies will be presented in a separate paper ^£ 
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